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1. Introduction 



Non-Gaussianity [Q has become a very important probe to the physics in the early universe. It is 
helpful to figure out the mechanism for generating primordial curvature perturbation. A well-defined 
non-Gaussianity is the so-called local form non-Gaussianity which says that the curvature perturbation 
can be expanded to the non-linear orders at the same spatial point 



3 » ^ 9 



C(X) = C 9 (X) + -fNL(g(x) + ^9NLQ{*) + 



where Cg is the Gaussian part of curvature perturbation, /jvl and g^L are the non-Gaussianity param- 
eters which characterize the sizes of local form bispectrum and trispectrum respectively. Single-field 
inflation model predicts f NL ~ 0(n s - 1) ~ 0(1(T 2 ) @. A convincing detection of a large local form 
non-Gaussianity is the smoking gun for the multi-field inflation in which the quantum fluctuations of 
the isocurvature field contribute to the final curvature perturbation on the super-horizon scales. See, 
for example, @ |, § §, |, § |, 0, 0, 0, ^ |l|, ^ [l§ 0, [l|, |l|, || etc. 

In the literatures the non-Gaussianity parameters are taken as a constant. However, recently 
there are some hints at a possible scale dependence of /nl come from the numerous observations of 
massive high-redshift clusters |2|, ||] which seems in greater abundances than expected from 
a Gaussian statistics [25, ^6f[ . In fact, may also be scale dependent. The scale dependences of 
/jvz, and qnl are measured by their spectral indices n^L and n 
which are respectively defined by 



9NL and their runnings a f NL and a gNL 



fNh{k) = fNL(k p ) 



k \ n lNL + 2 a fNL lD fc^ 



gNL(k) = 9NL{k p ) — 



\ n SNL + 1 a 9NL ln 7^ 



(1.2) 
(1.3) 



where k p is a pivot scale. The authors in [27] showed that Planck [28] and CMBPol [29] are able to 
provide a 1-a uncertainty on the spectral index of Jnl as follows 



An fNL ~ 0.1 



50 



fNL 



for Planck, 



ky 



:i-4) 
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and 



An 



Inl 



0.05 



50 



fNL y/f. 



for CMBPol, 



(1.5) 



sky 



where f s k y is the sky fraction. The effects on the halo bias from the scale dependent Jml are discussed 
in |3l| |. The studies on fingerprints of the scale-dependent g^i in CMB and large-scale structure 
are called for in the near future. 

Actually the scale-independent non-Gaussianity parameters are not generic predictions of inflation 
model. The large scale dependences of non-Gaussianity parameters can be obtained in the axion re- 
flation with different decay constants for different axion fields [32] and in the model with self- interacting 



canonical isocurvature field [33, 34 



37, 38, 39]. It is interesting for us to ask how the results 



should be modified for the non-canonical isocurvature field. In particular, we are curious that the 
scale dependence is enhanced or suppressed by the sound speed. For simplicity, we focus on the case 
with a DBI isocurvature field. 

Our paper is organized as follows. In Sec. 2 we will discuss the dynamics and quantum fluctuation 
of a DBI isocurvature field. In Sec. 3 the spectral indices of /jvx an d gNL from a DBI isocurvature 
field will be calculated. Some discussions are contained in Sec. 4. 



2. Dynamics and quantum fluctuation of a DBI isocurvature field 

Brane inflation [^] is considered to be a popular inflation model embedded into string theory. A more 
realistic setup is that a D3-brane moves in the internal six-dimensional Calabi-Yau manifold which 



contains one or more throats in Type IIB string theory [41]. For simplicity, we consider a D3-brane 
which is mobile in a throat whose metric is given by 



ds 2 = hT^g^dx^dx" + h 1 ' 2 [dr 2 + b 2 (r)d0 2 + 



(2.1) 



Here 9 is an angular coordinate which is transverse to the radius direction r and b(r) is the radius of 
throat at r. In this paper, we choose signature (—,+,+, +). The action for the mobile D3-brane is 
given by 



S 



J d 4 x yf-det [h-^ig^ + hd^rd v r + hb 2 d^9d v e)] - h~ l - J d 4 x^V(r,9). (2.2) 



It is convenient to define two canonical fields in the slow-rolling limit as follows 



T 3 r, 
%b9. 



(2.3) 
(2.4) 



Here we consider a 4>, and then (j) and a are taken to be the adiabatic and entropic direction during 
inflation respectively. 

In the limit of f(j) 2 <C 1 and fa 2 <C 1, these two-field inflation is reduced to the canonical case, 
where f = h/T^. In this paper we focus on another limit in which 



(2.5) 



and 



a 2 IS 2 « c 2 . 



(2.6) 
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From the above action the equations of motion for (f> and a are roughly given by 



3H(1 - K/3)<j) 



2f 



1 df dV(<t>,(j) 



3Ha ~ c, | 3 — ii a 



where 



Vb 



Hc s 

_b_ 

Hb 1 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



and both k and rji, are assumed to be much smaller than unity. For simplicity, the cross-coupling 
between cf> and a is assumed to be negligibly small, and Eq. fl2.8|) is valid when 



V"{a) 



3H 2 



< c 



-i 



(2.11) 



This is the slow-roll condition for the isocurvature field a. Now the dynamics of a becomes 

3H& ~ -c s [mV + V'(a)] , (2.12) 

where 

m 2 = -3i] b H 2 /c s . (2.13) 

The slow variation of the radius of throat induces an effective mass for the isocurvature field along 
the angular direction. 

The quantum fluctuations of <f> and a have been well studied in [42]. See also [18, 43, 44]. Here we 
don't want to repeat the computations. We will only briefly recall the results in [^]. The canonically 
normalized quantum fluctuations of (f) and a are respectively given by 



V<f, 



3/2 
Cs 



5(f), v a 



z5a, 



(2.14) 



whose Fourier modes corresponding to the Minkowski-like vacuum on very small scales are given by 



' e~ ikc °x [ 1 



where 



y/2kc~. 



X 



kc s x 



dt 



(2.15) 



(2.16) 



is the conformal time. The perturbation mode of k exits horizon at the time of c s k = aH. Therefore 
the power spectra for 8<f> and 5a are 



2vr J 
H 

27TC.O 



(2.17) 
(2.18) 
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The amplitude of 8a is amplified by a factor of l/c s compared to the canonical one. 

Before closing this section, we want to estimate the typical value of a during inflation. The 
quantity of (a 2 ) coming from its quantum fluctuations is 



1 



(2tt)2 



\Sa k \ 2 d 3 k 



1 



27T) 



1 



+ 



H 2 



2a 2 k 2k 3 c 2 



d 3 k. 



(2.19) 



Considering that the physical momentum p is related to k by p = k/a = e k, we have 



<- 2 > 



d 3 p fl H 2 



(2vr) 3 J p \2 2p 2 c 2 ) ' 



(2.20) 



The first term is contributed from vacuum fluctuations in Minkowski space and it can be eliminated 
by renormalization. In addition, the physical momenta for the modes of quantum fluctuations we 
concern in the time interval between and t are those from c7 l H to c7 1 He~ Ht . Therefore 



<- 2 > 



H 2 



c7 L H 



dp H 3 



^ 2 C 2 J c ^He-Ht p 4:TT 2 C 2 



(2.21) 



It indicates that the quantum fluctuation of the field a in the inflation epoch can be modeled by a 
random walk with step size H/2ttc s per Hubble time. However the vacuum expectation value of a 
cannot go like t for t — > oo if a has a potential. For example, let's consider the potential of a is 



V(a)= l - m la 2 + Xa n . 



From Eq. (2.12), the dynamics of a is described by 

3H& = -c s m 2 a(l + n\a n ~ 2 /m 2 ), 

where 



2 2,-2 



(2.22) 



(2.23) 



(2.24) 



Similar to 45], combining the contributions from quantum fluctuation and classical equation of motion, 
we have 



d{<7 2 



H 6 



2m" 



dt 4vr 2 c? 3H 



1 + ^ (a 3)»/2-l 



(2.25) 



The above differential equation approaches a constant equilibrium value, namely 

3H 4 1 



87r 2 m 2 c3l + §s' 



where 



s = 2A(cr ) ' /m . 

The typical value of a during inflation is a* = (a 2 ) ~ c s 3 ^ 2 which is enhanced for c s <l. 



(2.26) 



(2.27) 
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3. Scale dependence of local form non-Gaussianity parameters 



In this paper we consider that the curvature perturbation is generated by the isocurvature field a at 
the end of inflation or deep in the radiation dominated era, and the curvature perturbation can be 
expanded to the non-linear orders by using the 5N formalism 



C(*/,x) = N a (t f ,ti)5a(ti,x.) + ^N aa 5a 2 (ti,x) + -iV^&r 3 ^, x) + -, (3.1) 

where N a , N tfTa and iV )OW are the first, second and third order derivatives of the number of e- folds 
with respect to a respectively. Here tj denotes a final uniform energy density hypersurface and tj 



labels any spatially flat hypersurface after the horizon exit of a given mode. Similar to |36|, 38], ti is 
set to be t*{k) which is determined by k = aft*)!!* for a given mode with comoving wavenumber k. 
H* denotes the Hubble parameter during inflation from now on. 

From Eq. ( |3.lD , the amplitude of curvature perturbation is given by 

The amplitude of Gravitational waves perturbation only depends on the energy scale of inflation as 
follows 

A| = ^ (3-3) 

Here we work on the unit of M p = 1. The scale dependence of gravitational wave perturbation is 
measured by tit which is defined by 

n ^d^k=- 2 ^ (3 - 4) 

where 

For convenience, the tensor-scalar ratio tt is introduced to measure the amplitude of gravitational 
waves: 



Comparing Eq. (|3.1| ) to ([L . 1|) , the non-Gaussianity parameters are given by 



and 



= 6lvW ( } 



25 N aa(7 (U) 

9NL = mW ( } 



Following the method in 36, pq[ , we introduce a new time t r {> t*) which is chosen as a time soon 
after all the modes of interest exit the horizon during inflation and keep it fixed. The value of a at t r 
is related to that at time t* and the time t* by 

d(J - ^ C ^-dt, (3.9) 



y c '» J u 3H(t) 
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where 



Ves = -rh 2 a 2 + V(a). 
The equation of motion in Eq. ( |2.12| ) can be rewritten by 

3H& = -c s V^ s . 



Therefore we have 



da r 
da* 
da r 
dt* 



(3.10) 

(3.11) 

(3.12) 
(3.13) 

(3.14) 

which implies that F(a r ) is scale independent. 1 Taking into account that oy is a function of a*, we 
have 

da r dN(a r ) 



K ff K) 

C s {t*)V^{a r 



Considering 



and Eq. ( |3.11| ), one finds 



—F(a ) = (cr d°V + 9a r , 

dt* da r da* dt* ' 



d 



-F(a r ) 



d 



-F(a r ) = 0, 



da* da r ' 



d 2 a r dN(a r ) ( da r \ 2 d 2 N(a 
da 2 da r V da* J da 2 



N f t ) - d 3 a r dN(a r ) | ^da r d 2 a r d 2 N(a r ) | /da r \ 3 d 3 N(a r 
dal da r da* da* da 2 \da*J da 3 . 



(3.16) 
(3.17) 
(3.18) 



Since ffla ^1 and 

Gov ' <3(T^ 



are scale independent, one obtains 

dlnN )CT (t* 



dink 
dlnN, aa (t*) 

din A; 

d\nN^ aa (t*) 



c s r] a 



2c s 7] aa + c s rj 3 



N a (t*) 



3c n i 3c m jW^j | c * ^ 



din A; " " N i<7(7(7 (U) 

where the slow-roll equation of motion for cr is considered and 



>lo 



V"z(°*) _ _ K'ffV) l - V 



3 ^2 , 3ff a > v. 3ff 2 • 

From the above results, the spectral indices of A^, /tvl and g^L are respectively given by 



,(4), 
'off 



n s = 1 + 



din AS 



din A; 



1 - 2ch -2k + 2c s rj a 



(3.19) 
(3.20) 

(3.21) 
(3.22) 
(3.23) 



1 More precisely, the perturbation mode with k exits horizon when c s k = a*H*, and then dlnfc = (1 — en — K)H,dt* 
H«dt«. 
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din I/. 



NL 



dink 



c s r]3 



(3.24) 



and 



dln\g NL \ 



cHnA; ~~ S,J iV CT(JCT (i*) N i<Taa {t if/ 
We see that the spectral indices of both f^L and are suppressed by a factor of c s for c s <C 1 
Similar to [36, pg|j , the running of spectral indices of /jvx and g^L are defined by 



(3.25) 



<iln A; 



n 



9N L 



9NL 



dink 



(k + 2e# - c s r/ CT0 
(k + 2e H - 1c s r] aa )n 



n 2 



9nl 



n 



9NL 



1 



where 



9NL 



2/^i(c s r/ <TCT - c s r/ 4 + n /jVi )n /js 



3H 2 V" r 



3^yW • 



(3.26) 



(3.27) 



(3.28) 



The spectral indices of /jvx and are nice parameters to characterize the scale dependences of 
these two non-Gaussianity parameters only when nf NL and n gNL are much less than unity. 

3.1 DBI isocurvature field with polynomial potential 

In this subsection we will consider a simple example in which the DBI isocurvature field has an effective 
polynomial potential 

1 



V eS (a) = -m 2 a 2 + \a n , 



and then 



m 

3H 2 



1 + n(n — 1)- 



r] aa n(n - l)(re - 2)s/2 
cr* 1 + n(n — l)s/2 
Combing the normalization of curvature perturbation, we find 

, . fAT . 5 ^H*/2TTc s n(n-l){n-2)s/2 



6An 



(7* 



1 + n(n - l)s/2 



(3.29) 

(3.30) 
(3.31) 

(3.32) 



Here An is normalized to be 4.96 x 10 5 by WMAP in 47]. The slow-roll condition for a is c s r] aa <C 1. 
H*/2tvc s is the amplitude of quantum fluctuation of a and it should be less than <r*. For c s "(] aa ~ 10 -2 , 
H * l ^ Cs ~ 10" 1 and s > 10 _1 , we have K^/tvlI ~ O(10) which is detectable by PLANCK. In this 
case, if c s < 10~ 2 , r\ aa > 1 which implies that the effective mass of a is not less than the Hubble 
parameter during inflation. 

Taking into account the typical value of a in Sec. 2, the above equation becomes 

nj NL fNL = sign(N >a )(c s r] aa ) 3/2 d(s), (3.33) 

where 

5\/2 / 



d(s) 



l + ns/2 n(n-l)(n-2)s/2 



(3.34) 



6A n y 1 +n(n- l)s/2 l + n(n-l)s/2 

which is illustrated in Fig. |l]. For n = 4 and s>l, d(s) ~ 2.7 x 10 4 , the scale dependence of /jvx can 
be detected by PLANCK if c s r] aa > 0.003. 
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Figure 1: The function of d(s). The solid and dashed curves correspond to n — 4 and n — 8 respectively. 
4. Discussions 

In this paper we calculate the spectral indices of f^L and g^L and their runnings generated by a DBI 
isocurvature field on the super-horizon scales. We find that the indices are suppressed by the speed of 
sound c s . This suppression effect is not surprised and one can understand it from the effect of Lorentz 
boost. Here l/c s is nothing but the Lorentz boost factor. The time interval in the frame where brane 
is moving is dilated by a factor l/c s compared to that in the frame where brane is at rest. Because the 
coordinate a is transverse to the motion direction of brane, the value of a does not change. Therefore 
the velocity of a is suppressed by a factor c s . That is why there is a factor c s on the right hand 
side of Eq. (|2.12| ). On the other hand, the scale dependences of the non-Gaussianity parameters come 
from the small variation of isocurvature field a from tk to tk+dk, where tk and tfc+dfc correspond to the 
time when perturbation modes k and k + dk exit horizon during inflation respectively. Therefore the 
variation of isocurvature field a is suppressed by c s , and hence the spectral indices of /jvl and g^i 
are suppressed by c s as well. 

Even though the scale dependences of non-Gaussianity parameters generated by the DBI isocur- 
vature field are probably detectable, its effective of mass is required to be comparable or larger than 
the Hubble parameter during inflation if c s -C 1. 

There are so many models which can produce a large local form non-Gaussianity and an accurate 
measurement of Jnl cannot help us to distinguish them. But once the scale dependences of the non- 
Gaussianity parameters are detected, we can get more information about the isocurvature field, such 
as how it interacts with itself. 
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